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We study the noise correlations of one-dimensional binary Bose mixtures, as a probe of their quan- 
tum phases. In previous work [23], we found a rich structure of many-body phases in such mixtures, 
such as paired and counterflow superfluidity. Here we investigate the signature of these phases in 
the noise correlations of the atomic cloud after time-of-flight expansion, using both Luttinger liquid 
theory and the time-evolving block decimation (TEBD) method. We find that paired and coun- 
terflow superfluidity exhibit distinctive features in the noise spectra. We treat both extended and 
inhomogeneous systems, and our numerical work shows that the essential physics of the extended 
systems is present in the trapped-atom systems of current experimental interest. For paired and 
counterflow superfluid phases, we suggest methods for extracting Luttinger parameters from noise 
correlation spectroscopy. 



I. INTRODUCTION 

In recent years, the study of noise correlations has been 
established as a way of probing ultracold atom systems 
P~H5]- First proposed in Ref. [3], noise correlation spec- 
troscopy has been discussed as a way of measuring cor- 
relation functions of cold atom systems, such as pairing 
or density order [5HS]- In the experiments reported in 
Refs. [IHS], a cold atomic gas is first held in a trap, and 
then released from it by turning off the trapping poten- 
tial. The noise correlations are measured as the spatial 
correlations of the density in the fully expanded atomic 
cloud. If atomic interactions during the expansion can 
be ignored, one can use the noise correlation measure- 
ments to infer momentum space correlations in the ini- 
tial state. Such an analysis has been used to demonstrate 
the phase transition between superfluid (SF) and Mott 
insulator (MI) states [U [3] , as well as the formation of 
fermionic pairs [5]. 

The experimental realization of quasi-one dimensional 
many-body systems with ultra-cold atoms in optical lat- 
tices has been reported in Refs. |10H18| . Characteristic 
features of such systems include fluctuating and compet- 
ing orders. In contrast to higher dimensional systems 
which exhibit long-range orders, ID systems typically 
display only quasi-orders, that are characterized by the 
algebraic decay of the correlation function of the order 
parameter. In Refs. [BJIS], noise correlations were shown 
to be an effective probe of such orders in ID Fermi sys- 
tems, for both one- and two-component systems. Similar 
studies have been done for ID bosonic systems, either 
in the hard-core limit |19| or using Luttinger liquid (LL) 
theory [5 . In Ref. [5] , the signature of condensates and 
quasi-condensates was discussed in detail. 

Noise correlations can also be used to study the phases 
of binary bosonic mixtures. In such mixtures, two addi- 
tional orders beyond SF and MI are potentially present, 
first studied in Ref. [5D]. If the inter-species interaction 
is attractive, bosons of different species can form a paired 
superfluid (PSF) state. If the interaction is repulsive 
and the system is confined in a lattice at half-filling, the 
bosons can form particle-hole pairs, called "anti-pairs". 



Such anti-pairs can then form a counter-flow superfluid 
(CFSF) state. In addition, a charge density wave or- 
der (CDW) can coexist with the three superfluid orders, 
which is the defining feature of supersolidity. Numerous 
examples of such order have been given in Refs. [211 122| . 
In Ref. |2S] we established the phase diagram of a binary 
mixture exhibiting SF, PSF, CFSF and MI orders, and 
we showed that each of the superfluid orders can coexist 
with the CDW order. 

We also showed in [23] that because the PSF and CFSF 
orders are the result of inter-species pairing, they do not 
provide a signature in the momentum distributions of the 
individual atomic species. In this paper, we show that 
noise correlation measurements provide distinctive sig- 
nals of both the PSF and CFSF orders. Ref. [24] shows 
that noise correlations characteristic of the PSF/CFSF 
orders can be observed even in a system of only four 
atoms. Here we calculate the noise correlation spectra 
from first principles, using the time-evolving block deci- 
mation (TEBD) method [55], which is supported by an- 
alytical calculations based on LL theory. We make ap- 
propriate comparisons between results for homogeneous 
and trapped systems. 

To evaluate the noise correlations, we first assume bal- 
listic expansion and long expansion time and define the 
noise correlations as the density correlations in momen- 
tum space, 

Gaa'(k,k') = (n a ,kn a ',k') - {n a ,k){n a >,k') (1) 

where a, a' are species indices (a, a' = 1,2), k, k' are 
momenta, and n a ,k an d n a / t k' are the occupation opera- 
tors in momentum space. We also consider the derived 
quantities C aa i(q) and D aa >(q), defined as 

Caa'(2 q ) = f dk ^ + t" { n~ 9 \ > ( 2 ) 



and 



D aa , (2q) 



dk 



[n a (k + q))(n a >(k - q)) ' 

(n a (k + q)n a >(q - k)) 
(n a (k + q)){n a ,{q - k)) ' 



(3) 



Each of these quantities can capture the main features 
of the noise correlations for particular types of order and 



2 



can be directly measured in experiments [3 J. We will 
present all our results first in the form of Q aa '{k, k') and 
then use C aa '(q) and D aa i(q) to highlight the key fea- 
tures. 

This paper is organized as follows: in Sec. [TTJ we ex- 
plain our model and the different quasi-orders present in 
it; in Sec. |III| we show our LL results and predict the 
generic signature of the noise correlations for different 
orders. In Sec. |IV| we present our numerical calculation 
of noise correlations for both homogeneous and trapped 
systems. In Sec. [V]we conclude. 



II. NOISE CORRELATIONS AND 
QUASI-ORDERS 

We work in an approximation in which binary Bose 
mixtures in optical lattices are described by a two- 
component Bose-Hubbard model [23, 26J. The Hamilto- 
nian for M atoms of each species confined to an optical 
lattice with N sites is given by 

N-l N 
et=l,2 j=l j=l 



N 



o=l,2 j = l 



(4) 



We denote the different atomic species with the index a — 
1,2, and the lattice site with index j. We assume that 
the two species have the same average filling factor, v — 
M/N < 1, the same intra-species interaction U > and 
hopping parameter t > 0. The inter-species interaction is 
given by U\2- The operators b\ ^ and b a j are the creation 
and annihilation operators for atoms of type a and site i 
and n a j = b^ a jb a j are the number operators. 

In Ref. |22], we found that there are four different 
regimes in the phase diagram (besides the collapsed or 
phase-separated regime that occurs at large | CT - 12 1 ) : the 
superfluid (SF), the paired superfluid (PSF), the coun- 
terflow superfluid (CFSF) and the Mott insulator (MI) 
state. In addition, each of the superfluid orders can co- 
exist with a charge density wave (CDW) order. The ex- 
istence of any such order is determined from the asymp- 
totic behavior of the correlation functions of the cor- 
responding order parameters. Specifically, the single- 
species superfluid (SF) has the order parameter b a (x) 
(a = 1,2) and its corresponding correlation function 
G(x) = (bl(x)b a (0)}; the paired superfluid (PSF) has 
the order parameter bx (35)62 (&) and its corresponding 
correlation function Rs(x) — (b\(x)b\(x)bi(0)b 2 {0)) ; the 
counter-flow superfluid (CFSF) has the order parame- 
ter b\(x)b2{x) and its corresponding correlation function 

R A {x) = {b\tx)b 2 {x)bi (0)4(0))- The CDW order param- 
eter is n a (x) (a = 1,2) and the corresponding correlation 
function R na {x) = (n a (x)n a (0)) . The asymptotic be- 
havior of the correlation functions at large x is listed in 
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Table I: Definitions of MI, SF, CFSF and PSF orders in terms 
of the large x behavior of the correlation functions Rs(x), 
Ra(x), and G(x), R n , a ( x )- A: algebraic decay of the form 
x~ a ; E: exponential decay of the form e~^ x . A correlation 
function is said to exhibit quasi-order when it is subject to 
algebraic decay with a < 2. In this system, the algebraic 
decay for Rs, Ra and G always has a < 2, while R n ,a can have 
a > 2. CDW quasi-order exists only when R„ t a is described 
by a < 2. 



Table. |T]for the different phases. 

We calculate the noise correlations from the four-point 
correlation function: 

N 

g aa ,(k,k') = £aa>{h,h,h,h)<? [kj " +k ' hi] 

J'li,i2,i3,j4=l 

-{n a {k)){n a ,{k')), (5) 

where j 12 = ji - j 2 , j 3 4 = h - H and ^ is the four- 
point correlation function, 



Caa'{h,h,h,H) = { h a, n Kh b l',j 3 h a'u)- (6) 

It is easy to see that the correlation functions Rs, Ra 
and i?„ j0 are the special cases of C aa >, 



Cl2(jl,j2,j2,jl) = RA(jlj2), (7) 

The noise correlation Q12, therefore, contains the 
Fourier transform of Rs and Ra, 



and 



as 



3A 



■h) 



(k,k')=J2 R A(jl,32) 



,i(k-k')(j 1 -j 2 ) 



(8) 



(9) 



and Q aa contains the Fourier transform of R n . a , 



32) 



(10) 



If RsUuh) decays as \ji - j 2 \ 1/Ks , we find that 
g s scales as \k + k'\~ x / Ks . Similarly, if Ra decays as 
b'l - J2\~ 1/Ka , 3 a scales as |fc - k'\-^ KA . For the PSF 
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phase, we find that g s (k,k') is the dominant term of 
Q\ 2 {k,k') with a strong peak around k = —k'. This 
peak is the signal of the PSF order. Similarly, for the 
CFSF phase, we find that the function gA(k, k') becomes 
dominant around k = k! in Qi 2 (k, k'). The peak around 
k = k' is the signal of the CFSF order. These remarks are 
made to give the reader an intuitive interpretation of the 
relationship between the noise correlations and the long- 
range orders. In the following section, we will explain 
the calculation of the noise correlations via LL theory 
and show that the features mentioned above are indeed 
reflected in the LL calculation results. 



The parameters K$,a are the Luttinger parameters of 
the symmetric/anti-symmetric sector. To calculate the 
noise correlations away from the SF regime, we take the 
limits Ks,a — > to describe the phases in which either 
or both Rs/a have short-ranged correlations (exponen- 
tial decay) . This approximation corresponds to the limit 
that the length scale of the exponential decay is much 
smaller than any other length scale of the system. Our 
calculation could be extended in a straightforward way 
to include a finite decay length of the exponential decay. 

We start out by calculating {n a j~) for small momentum 
k w 0, for which the Bose operators are given by b a ~ 



For (n a k) we find: 



III. LUTTINGER LIQUID APPROACH 

In this section we determine the generic behavior of the 
noise correlations using a Luttinger liquid approach. This 
formalism has been applied to one-dimensional Fermi 
systems in Ref. [B], and additionally to single-species 
bosonic systems in Ref. |5_ , where a detailed description 
of these calculations was given. Here, we use an anal- 
ogous derivation for the case of a bosonic mixture. We 
outline key steps of the derivation, but refer the reader 
to Ref. (5] for a detailed description of the method. 

As described in Ref. [33] , we switch to a continuum de- 
scription, in which the single particle operators are given 
by b a (x). We then use a bosonization identity |271 E5] 



b a (x) = [n + n a (x)} 1/2 Y, e2 " nea(x)e ^ Jx) : 



(11) 



where n is the real space density, related to the filling 
factor by n — v/cll, where a^is the lattice constant, and 
m is an integer summation index. For future reference, 
we note that the Fermi wavevector of an equivalent sys- 
tem of fermions, fc^, is given by kp — nn. Although this 
paper describes a bosonic system, we find that the Fermi 
wavevector occurs naturally in a number of contexts. For 
example, O a is given by Q a (x) = kpx + 9 a {x), where 
9a( x ) — 71 J dy\l a {y). n a describes the low-momentum 
density fluctuations of species a and 4> a (x) is the phase 
field of species a. 

We calculate the noise correlations at the Gaussian 
fixed point, corresponding to the SF phase. Here, the 
system separates into symmetric and anti-symmetric de- 
grees of freedom, defined as 9s. a — (#i ± &2)/V2 and 
4*s,A — (4>i ± </ ) 2)/v / 2- The action can be written either 
in terms of the phase fields 



S = 



E 

j=S,A 



d 2 , 



K 
2n 



3 -[{d V] r^) 2 + {d x ^) 2 ] , (12) 



or in terms of the fields 



'S,A 



E 

j=S,A 



d 2 r, 



2nK, 



\{d VjT eif + (d x 3 ) 2 } .(13) 



The velocities v$.a are the phonon velocities of the 
symmetric/anti-symmetric modes, and r$,A = (vs.Ai~,x). 



(n a ,k) 



dx 12 e^e~^^-^ 2 \ 



(14) 



where a (l) refers to 4> a (x\), and similarly for a (2), and 
X12 = x\ — X2- The correlation function ((</> a (2) — a (l)) 2 ) 
can be rewritten in terms of correlation functions for 
4>s.a- Using the Gaussian action above, we find 



(Ws M (2)-^ M (l)) 2 > = ^^log T i±^,(LV i 



2K s/A hS r 



where ro is a short-range cut-off. With that we find 

(n k ) ~n[dx 12 e ikx ™F(x 12 ), (16) 

where 

n*) =(^LY- (17) 



The exponent g is given by g — I/8X5 + 1/8Ka- Next 
we evaluate the expectation value (n^n^) along the same 
lines. We obtain: 

(ni, fc m,fc/ )~ n 2 j e ikxM ' x ™ F(x 12 )F(x 34 )A(18) 

where 

A — ( ( r l+ x ll){ r l+ x2 23) \ h (1Q\ 

and Eq. [18] is a volume integral over the spatial vari- 
ables x±2, X23, X34. The exponent h is given by h = 
— l/8Ks — 1/8Ka- We combine these expressions to get 
the correlation function Gn(k, k'): 



MM') 



n 2 I e 



ikxi2+ik ' X3i F(x 12 )F(x M )(A - 1). (20) 



For G\ 2 {k,k') we proceed analogously, and find h — 
— \/8Ks + 1/8Ka- For the finite-size systems that we 
treat here, we evaluate these integrals numerically, by 
choosing a finite length L of the system, and by replac- 
ing each spatial variable x by (L/2tt) sin(27ra;/L) (see Ref. 
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(a) MI (b) SF (c) PSF (d) CFSF 



Figure 1: Noise correlations in different phases, derived from Luttinger liquid theory. In the MI state (column (a)), 5-function 
like correlations along k = k' in Qi\{k, k') are visible, whereas Q12 nearly vanishes. In the SF state (column (b)), with Luttinger 
parameters Ka = 1-03 and Ks = 0.96, we find various contributions in Qn, especially S- function along k = k' . In Fig. [2] we 
show the contour plots for C/11 and C/12 for the same state, where we can see the negative correlations at k = and k' = 0, as 
well as pairing correlation along k — —k', which is similar to the single-species result in Ref. [5]. G12 shows similar features, 
but the bunching contribution is an algebraic peak, rather than a <5-function. In (c) we show an example for the PSF phase, 
with K A = 0.01 and K s ~ 1.3, in (d) an example for the CFSF phase, with Ks = 0.01 and K A ^ 1.2. In the PSF state, the 
inter-species correlation G\2{k, k') has strong correlations along k = — k', a reflection of pairing. In the CFSF state, the peak 
is formed along k — k' direction, an indication of the anti-pairing (particle-hole) formation in the CFSF state. 
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Figure 2: Noise correlations, Qn(k,k') (a) and Qi2(k,k') (b), 
in the SF state. The data used here are the same as those 
used in Fig. [I] (b) . We create the non-linear gray scales by 
plotting functions tanh(500<5n) (a) and tanh(10 5 C/i2) (b), in 
order to magnify the details around k = k' = 0. The labels 
of the color-bar reflect the values of Q\\(k, k') and Gi2(k, k'). 
In these plots, we can clearly see the negative correlations 
between the quasi-condensate (k — and k! — 0) and the 
higher momentum states at the quantum depletion, as well 
as the anti-pairing correlation along k = k' and the pairing 
correlation along k = —k' . 



|29|). To compare with the TEBD calculations of a homo- 
geneous system in the next section, we choose the values 
of the Luttinger parameters, Ka and Ks to those ob- 



tained from the TEBD calculations and they are listed 
in Table. HU 

In Fig. [l] (b) and [2] we show an example for the SF 
regime. In the upper panel of Fig. [I] (b), we show 
Gn(k,k'), in the lower panel Gi2{k,k'). The Luttinger 
parameters are Ka — 1-03 and Kg — 0.96. The ratio 
L/tq was chosen as L/r = 20, corresponding to the par- 
ticle number of each species in the numerical example. 
The shape of Qn{k,k') is the same as the noise corre- 
lation function for a single bosonic SF, which was dis- 
cussed in Ref. |5]. It has the characteristic features of 
a superfluid: positive correlations along k = — k' , which 
indicates pairing correlations; negative correlations for 
the axes k = and k' — 0, indicating the negative cor- 
relations between the quasi-condensate and the higher 
momenta due to pair fluctuations; and bunching correla- 
tions along k — k! . For Qvi{k, k') we find qualitatively a 
similar shape, with the main difference, that the bunch- 
ing along k — k' does not have a ^-function contribution, 
but only algebraic terms. We note that for a system of 
two non-interacting species, i.e. for U12 = 0, Ks = Ka, 
and Q12 vanishes. 

In Figs, [l] (c) and (d) we show the noise correlations 
for the PSF and the CFSF phase, respectively. For the 
PSF example, the Luttinger parameters are Ks = 1.3 
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and Ka = 0.01. For L/tq we again pick L/tq = 20. 
For the CFSF phase, the parameters are Ka = 1.2 and 
ifj = 0.01. In the PSF regime, we find a strong pairing 
signature in Q12, similar to the pairing signature in Fermi 
mixtures [51 |B]. In the CFSF example, an strong anti- 
pairing signature is found in Qi2- 

We can obtain the functional form of these signatures 
in the limit L — > 00, by applying similar arguments to 
what has been given in Ref. [S]. In the PSF region, 
We rewrite the noise correlation integral in terms of z = 
(xi2 - x 34 )/2, h + = (X14 + £ 2 3)/2 and h_ = (x u - 
cc 2 3)/2. We then note that for Q 12 and for Ka — > 0, the 
exponent h = —l/8Ks + 1/8Ka diverges. This enforces 
the integrand to be negligible away from z, h + sa 0. Thus 
the integral evaluates to 



\k + k'\- 1 ' K s 



(21) 



This is the shape that would be approached in an infi- 
nite system by the noise correlations shown in Fig. [l] 
(c), lower panel. The deviation from the pure power 
law is due to the finite size of the system. With simi- 
lar arguments one can show that in the CFSF regime the 
inter-species noise correlation approaches 



\k-k'\~ 1 ' KA , 



(22) 



for L — > 00. Again, the deviation from a pure power law 
is due to the finite size of the system. Furthermore, one 
can show that for both PSF and CFSF orders, Gu(k,k') 
approaches S(k — k'), in the limit of infinite size. Equa- 
tions [21] and [22] show that there is a simple relationship 
between Q12 and K$ in the PSF regime and Q12 and Ka 
in the CFSF regime. This suggests that a careful mea- 
surement of Q12 can be used to extract the value of the 
Luttinger parameters appropriate to the system. This 
is further confirmed by our numerical calculations for a 
trapped system, where we show that the algebraic rela- 
tionship described by Eqs. [2T] and [22] remains valid in 
the presence of a harmonic trap. We discuss prospects 
for experimental determination of Luttinger parameters 
in Sec. ITTCl 

The MI result in Fig. [T] (a) is obtained by setting 
both Ka and K$ to 0.01. In this case, the ground state 
closely approximates a simple product of MI states of 
each species. Thus, Qu approaches a (^-function, whereas 
Qi 2 nearly vanishes. 

Next we calculate the noise correlations for the case 
k w and k! w 2k f, where kp is the Fermi wave vector 
defined above. Essentially the same calculation can be 
done for k! ~ —2kp, and k! sa and k sa ±2kp. n a ^ 
is still given by the expression (16 1, but n a .fc< now needs 



to be calculated with the operator representation b{x) 
y/nex-p(2iO(x)) ex-p(i<f)(x)). With that we find 



(n q '+2k P ) ~n / dx^v^T'ix^), 



(23) 



where J r '(xi2) has the same form as before but with an 
exponent g' = l/8K s + l/8K A + {K s + K A )/2. The noise 
correlations take the form 



Qu(k,q') 



n 2 I e ikx 12 +iq'x 34 



T{x l2 )T'{xu){A-l). (24) 



We therefore note that around the points k sa and 
k! sa ±2kp, and k sa ±2kp and k' sa the integrand is 
multiplied by a contribution that is of the form of the 
integrand of the static structure factor 



S(q) 



(K s +K A )/2 



(25) 



which can create cusps in the noise correlation when the 
system is in the CDW regime. These cusps are found in 
our numerical calculations and are discussed in the next 
section. 



IV. NUMERICAL CALCULATION 

The calculation of noise correlations is based on the 
ground state generated by the time-evolving block deci- 
mation method (TEBD). This method has been used to 
generate the ground state of many ID models [5S]. In 
this method, the Hilbert space H is decomposed as 



H 



(26) 



Here, I refers to the /th lattice site, N is the number of 
sites, and H; is the local Hilbert space at site I with local 
dimension d, independent of I. Any state in H is 
represented as 



I*) 



E 



\h)\h)-\3M), (27) 



where 



X X 



y y ... v A [1] r [llji A [2] r [21j2 a [3] ••• 

ai — 1 OL2 — 1 qm-i — 1 

x [M-l] v [M-l]j M -i\[M] r [M]j M \[M+l]f2H) 
A/1 a M -i om-ium "a M a M -ia M a M + l \ > 

The variables Aq, and xi are the Schmidt coefficients and 
rank of the Schmidt decomposition of l^) at site / and 
is a rank-three tensor. Further detail on this method 
is provided in the appendix of the previous publication 
[23 . Here, we limit ourselves to stating values of param- 
eters and particular methods of calculation. In this work, 
we set the Schmidt rank \ = 100 and the local dimension 
d = 5. We use imaginary-time propagation to generate 
the ground state. After obtaining the ground state, we 
calculate the correlation functions, Ra, Rs, R n .a and 
G, and determine the quasi-long range order present in 
the system based on the relationship shown in Table. [I] 
Furthermore, we can extract the value of the Luttinger 
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Parameter setting 


Order 


Luttinger Parameters 


(a) U12/U = 0.01, v = 1 


Mott Insulator (MI) 


Ka = Ks = 


(b) U12/U = 0.01, v = 0.5 


Superfluid (SF) 


K A ~ 1.03, K s ~ 0.96 


(c) U 12 /U = -0.11, v = 0.5 


Paired Superfluid (PSF) 


K A = 0, K s ~ 1.3 


(d) U12/U = 0.11, v = 0.5 


Counterflow Superfluid with (CFSF) 


Ka ~ 1.2, K s = 


(e) U12/U = 0.26, v = 0.2 


Superfluid with charge density wave (SF/CDW) 


K A ^ 1.4, Ks ~ 0.57 



Table II: The parameters used in the numerical examples and the Luttinger parameters extracted from the algebraic fit of 
correlation functions, Ra and Rs- The Luttinger parameters are set to zero when the correlations decay exponentially. The 
hopping parameter t is 0.02(7 for all cases. The parameters are chosen to represent different orders that can exist in this system. 



parameters, Ka and Ks, from the numerically calculated 
correlation functions [55] • We use these parameters in a 
LL calculation to compare the numerical and the analyt- 
ical results. 

The main challenge of determining the noise correla- 
tion functions is the high computational cost of calcu- 
lating the four-point function, £a,a,'{j\i02iH-,ji) (Eq. |6j, 
which is estimated to scale as x 3 ^ 3 ^ 4 - For the system 
sizes used in this paper, we use parallel computing algo- 
rithms to speed up the calculation by parallelization the 
computation of C aa ' along the indices ji. 



A. Homogeneous system 

In this section we discuss the numerical results for noise 
correlations of a homogeneous system of 40 lattice sites, 
subject to the hard-wall or "open" boundary condition, 
in which the wave function is required to vanish on the 
fictitious sites of index and N + 1 implied by Eq. [4] 
We consider five parameter sets listed in Table [TTJ rep- 
resenting different regimes of the phase diagram of ID 
Bose mixtures. 



Superfluid and Mott insulator 

For the Hamiltonian of Eq. [4j in the non-interacting 
case, U12 = 0, SF and MI are the only two possible or- 
ders. In the interacting case, SF and MI orders are still 
encountered, when the inter-species interaction is weak. 
For the Hamiltonian of Eq. [4] with t -C U, the MI state 
exists for any \U\2\ ^ U, until the occurrence of collapse 
(U12 < — U) or phase separation (U12 > U). The SF state 
however exists only when \Ui2\ <C U . In either of SF or 
MI phases, the quasi-order is formed in each individual 
species and the cross-species correlation is weak. 

For the MI state (Fig. [3] (a)) we find that Qu(k',k) 
shows strong correlations along the direction k' = k, in 
agreement with the LL theory result shown in Fig. [T] 
(a). We also find that the correlations along k' = k are 
not uniform and that the peak along k' = k resembles 
a Lorentzian distribution in k imposed upon a constant. 
This Lorentzian is due to the characteristic scale of the 
correlation functions. This contribution was ignored in 



the before-mentioned approximation in the LL calcula- 
tion, but could be included in a straightforward man- 
ner. The cross-species noise correlation, Qn{k,k'), on 
the other hand, is essentially zero, indicating the absence 
of cross-species correlations in the MI state. For the SF 
state (Fig. [3j (b) and Fig. |J, we consider the case where 
there is weak repulsion between the two species (Table 
|ll](2)). The Luttinger parameters are Ka = 1.03 and 
Ks — 0.96, which were extracted from the correlation 
functions Rs and Ra by numerical fitting. From the up- 
per panel in[3](b), we see that Q\\ has the characteristic 
features of a quasi-condensate [S]: the positive correla- 
tions along k = —k', which indicate pairing; the nega- 
tive correlations between k = and finite k' , as well as 
between k' — and finite fc; and a (^-function like corre- 
lation along k — k' , corresponding to bosonic bunching. 
The lower panel in Fig. [3] (b) , we see Q12 , which shows 
similar features, except for the 5- function along k = k' , 
which is "softened" into a power-law divergence and a 
slight negative value at k = k' = 0. For a system of 
two non-interacting superfiuids, i.e. U12 — 0, we have 
K s = K A , and Q X2 = 0. 



Paired superfluid and counter-flow superfluid 

We now discuss the noise correlations of the PSF and 
CFSF states. The noise correlation Gi2(k,k') is particu- 
larly important for these two phases, because it can verify 
the existence of PSF and CFSF orders. Unlike SF and 
MI states, PSF and CFSF states are characterized by 
order parameters that contain both species and there- 
fore cannot be reflected in any single-species observables, 
such as the single-particle Green's function, G a (x) or the 
single-particle momentum distribution (23]. The noise 
correlation function Gi2{k,k') measures the correlations 
between the momentum occupancies of the two species, 
and thus provides a direct probe of these orders. We 
have shown in the previous section, that the peak along 
k = —k' in Gi2(k,k') indicates the PSF order and that 
along k = k! indicates the CFSF order. These features 
are verified in our numerical calculation of Q\2 from the 
ground state. 

In Fig. [3] (c), we show the noise correlations in the 
PSF state. The parameters are listed in (c) of Tab. |n] 
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(a) MI (b) SF (c) PSF (d) CFSF 



Figure 3: Noise correlations for a homogeneous system of 40 lattice sites, calculated with the TEBD method. The frames (a) 
- (d) correspond to the examples (a) - (d) in Table |H| In (a), we show the noise correlations of a MI state. In the plot of 
Qu(k, k'), there is a strong correlation along the direction k = k', whereas the noise correlation function Gi2(k,k') essentially 
vanishes. In (b), we show the noise correlations of a SF state. Here, we can see the peak around k = k' corresponding to the 
S- function bunching peak predicted by LL theory (see also Fig. [T](b)). For C/12, we find negative value at k — k' = 0, which is 
different from the LL result (Fig. [TJ(b)). Other structures predicted by LL theory can be seen in Fig. [4j where Q12 and are 
plotted in a non-linear color scale to magnify the structures around k = k' = 0. In (c) and (d), we show the noise correlations 
of the PSF and CFSF state, respectively. In the PSF state (c), the inter-species correlation Q\2(k,k') has strong correlations 
along k — —k', a consequence of pairing (see also Fig. [TJ (c)). In the CFSF state (d), the peak is formed along the direction 
k — k! , an indication of anti-pairing in the CFSF state (see also Fig. IT]). 



The existence of PSF order is - as usual - determined by 
the behavior of the Rs(x) and Ra{x). Ra(x) decays ex- 
ponentially and Rs(x) algebraically with Luttinger pa- 
rameter Ks — 1.3. For the noise correlation function 
Qvi{k, A/), we find that a peak is formed along k = —k', 
which is a consequence of the pairing correlations. In Fig. 
[3] (d) , we show our numerical results for the CFSF exam- 
ple (d) in Table |nj Based on the behavior of the Rs(x) 
and Ra(x) we verify that the system is in a CFSF state 
with Ka — 1.2, and an exponentially decaying R$(x). 
For Q12 we find that a peak is formed along the diagonal 
direction, as a result of correlations of anti-pairs (b\b\). 
These findings are consistent with the predictions of LL 
theory (see Fig. [Tj). We note that Qyi(k 1 k') is enhanced 
in magnitude in the PSF and the CFSF phase compared 
to the MI and the SF phase, with a strongly altered func- 
tional form. 



Charge density wave 

In certain parameter regimes of the phase diagram, 
charge density wave (CDW) order can coexist with each 
of the three superfluid orders, SF, PSF and CFSF. In 



Sect. |III| we use LL theory to show that CDW order 
can be reflected in the function Q\\ and that the behav- 
ior of Q\\ around k = k' ± 2kp resembles the structure 
factor S(k). The reason for the resemblance can be un- 
derstood in a simple way, by recalling the definition of 
the structure factor 



S ^ = n £ e ~ tt0Wl) - K>K»" (29) 
h J2 

As mentioned in Sect. [IT] the density correlation function 
is "contained" in the noise correlations and the term , 

N 

E (^ 1 6iJ.^ A 6iJi>e i[ * (,1 -* )+fc ' Cfa - Jl)] . 

is part of the full sum that needs to be taken for Q\\. 
This term can also be written as a function of the density 
operator, rij = b^bj , as 

N 

E {n n n ]2 )e-^ k '- k ^-^ + 8(k - k')(n J1 )e' t( - k '' k) ' 1 
31,32=1 
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Figure 4: Noise correlations, Gu(k,k') (a) and Q\i(k, k') (b), 
in the SF state of a homogeneous system. The values of 
Gn(k,k') and Gi2(k,k') are exactly the same as in Fig. ^ 
(b). We create non-linear gray scales by plotting tanh(10<5ii) 
and tanh(2005i2) in linear scales. The labels of the color-bar 
reflects the values of Qu(k,k') and Qi2(k,k'). The features 
around k — k' = are magnified as a result of the non-linear 
scale. In (a), we find the features predicted by LL calculations 
j2j (a)). In addition, we can see a weak correlation at around 
k = k' ± 2/cf , where &f = v x tv/cil = 0.57r/a_L. This is 
where a strong correlation (cusps) will develop if CDW order 
is present. This feature can also been s how n in LL calcula- 
tions at around k ~ and k' ~ 2kp (Eq. 241. In (b), we find 



that the structures along k = k' is similar with the ones in LL 
calculations, however, the structures along k — —k' is nega- 
tive, different from the LL predictions (see also Fig. The 
difference may be understood as a result of different boundary 
conditions used for the finite-size calculations: the numerical 
calculations use a "hard-wall" boundary condition, whereas 
the LL calculations assume a periodic boundary condition. 
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Figure 5: Left : The correlation function Cn (q), as defined in 
Eq. [2] right: the structure factor S(k) for a quasi-supersolid 
state. The parameters are given in example (5) of Table. [TTJ 
The Luttinger parameters are Ka ~ 1.4 and Ks ~ 0.57. The 
filling fraction is v = 0.2, hence the "Fermi wave vector" kF 
is 7r x 0.2. At momentum 2&f, both quantities develop cusps, 
indicating the presence of CDW order. 



This shows that Qu and S(k) (Eq. 29 1 have the same 
Fourier transform of the density correlation function. If 
S(k) develops cusps at ±2kp, where kp = ttv [23J, when 
CDW order is present, we expect Qu to have similar 
cusps at k = k' ±2kp. In Fig. [5] we show one example 
of a quasi-supersolid (SS) state [21 , where CDW order 
coexists with SF order. The parameters are listed in (e) 
of Table |Tlj In the plot, the correlation function Cu(q), 
an integration of Qu (k,k') along the direction k = k! (Eq. 



xlO ' 
G (k,k') 




G 12 (k,k') 
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Figure 6: Noise correlations in the trapped system. The sys- 
tem size is 80 sites and t/U = 0.02. In (a), the system is 
in the SF state. The particle number of each species is 30, 
the trap frequency 8 x 10" 5 C7 and U 12 /U = 0.01. In (b), 
the particle number of each species is 40, the trap frequency 
1 x 10~ 4 U and U12/U = -0.11. The system has both MI 
and PSF orders. The MI state forms a plateau at unit-filling 
at the center of the trap and the PSF is formed at the edge. 
The PSF state at the edge causes the small peak along the 
k — —k' direction, similar to the one in (c). However, this 
peak is at a much smaller amplitude than the one shown in 
(c), where the whole system is a PSF state. In (c), the parti- 
cle number of each species is 20, the trap frequency 1 x 10~ U 
and U12/U = —0.11. The whole system is in the PSF state. A 
strong pairing correlation is formed along k = —k! direction. 
In (d), the particle number of each species is 30, the trap fre- 
quency is 8 x 10" 5 [/ and U12/U = 0.2. The system has both 
CFSF and SF order. The CFSF order forms a plateau at 
half-filling at the center of the trap and the SF state towards 
the edges of the trap. The CFSF order causes a strong anti- 
pairing (particle-hole) correlation along k = ^'direction. At 
the same time, the SF order adds to the "dips" along k = 
and k' = 0. 



|2|, is compared with the structure factor S(k) of the same 
state. In both functions, we can see cusps appearing at 
±2k F . 



B. Noise correlations in the trapped system 

We now discuss how the different types of order are 
affected by the presence of a trapping potential. To sim- 
ulate the effect of a trap, we add a harmonic potential, 
^0 ~ jc) (^1,2 + n 2,j) to the Hubbard Hamiltonian in 
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Figure 7: Correlations 612(g) and D\2{q) for the states that 
are described in Fig. [7] In (a), we show the behavior of 612(9) 
(Eq. |2| in SF, MI, PSF and CFSF states. The strong anti- 
pairing (particle-hole) correlations in the CFSF state gives a 
strong signal around q — in 612(g). This strong signal is 
also unique to the CFSF state and therefore can be used to 
detect to the CFSF order. In (b), we show the behavior of 
D 12 (q) (Eq. [3} in SF, MI, PSF and CFSF states. The strong 
pairing correlation in the PSF state is the reason for the high 
peak around q — in 612(5). This suggests measuring 612(g) 
is a good way of detecting the PSF order. 
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Figure 8: Noise correlation 611(g) and structure factor S(k) 
in a PSF/CDW state. The system size is 80 sites and there 
are 20 particles of each species (y = 0.25). The trap frequency 
is Q. — 10~ 5 U , the hopping t = 0.02C7 and the inter-species 
interaction is U\2 = — Q.11J7. The density at the center of 
the trap is roughly 0.45 per site and the cusps are developed 
around ±0.9-71. The inhomogeneity of a trapped system means 
that the "Fermi wave vector" Uf is no longer nu, where v is 
the average filling of the system. Instead, fci? can be evaluated 
as Tvn ce „t er , where n center is the density at the center of the 
trap, 



Eq. |4j where j is the site index and j c is the index at the 
center of the system. We then use the TEBD method 
to calculate the ground state. We also increase the sys- 
tem size to 80 lattice sites, and choose the total number 
of particles and the trap frequency to ensure that the 
boundary effect is negligible. 

One interesting feature of a trapped system is that dif- 
ferent orders can coexist in the trap. A well-known exam- 
ple is the MI plateau at the center of the trap surrounded 
by a SF at the edge [3U]. For repulsive inter-species in- 
teraction, we find coexistence of a CFSF plateau with a 
SF at its edge and a MI plateau with PSF at the edges 
for attractive inter-species interactions [23]. Despite the 
potential complication of coexistence of orders, we find 
clear signals for the pairing correlations of the PSF phase 
and the anti-pairing correlations of the CFSF phase. 

In Fig. [6] we show the behavior of Gi2(k,k') in four 
different cases, where the orders at the center of the trap 
are SF, MI, PSF and CFSF respectively. We find that the 
general behavior of the noise correlation in a trap is very 
similar to its homogeneous counterpart. In Fig. |6](c) and 
(d), G12 shows clearly the feature of pairing correlations 
in the PSF state and the anti-pairing correlations in the 
CFSF state. In addition, we see some minor features 
attributed to the coexisting orders. In the case of CFSF 
in a trapped system, we can see the "dip" along k = 
and k' = because of the coexistence with the SF order. 
On the other hand, in the case of a MI in a trap, we can 
see pairing correlations as a result of the residual PSF 
state at the edges. This pairing signal is much smaller 
than when the whole system is in the PSF state. 

To show that the peaks along k = k' and k = — k' 
in G12 are detectable in experiments, we also calculate 
Ci2(q) (Eq. [2| and Di 2 (q) (Eq. J3| for the four states. 
In the correlation Ci2(q) (Fig. una)), a high peak at 
q = only appears in the case of the CFSF state. This 



peak corresponds to the peak in G12 along k = k' in the 
CFSF state and is a reflection of the anti-pair correlation 
in the CFSF state. Similarly, in Di2(q). the high peak 
at q = only appears in the PSF state, as a result of the 
strong pairing correlations in the PSF state. A similar 
measurement has been performed for fermionic mixtures 
to detect the pairing of fermions [2]. 



In addition to the PSF and CFSF order, we also look 
for the signal of CDW order in Gn(k, k') in the trapped 
system. In a trapped system, the CDW order is more 
difficult to establish especially in the PSF and SF states, 
because the varying local density makes the "Fermi wave 
vector" im a spatially varying quantity. However, we 
can still see weakened cusps forming at the momentum 
roughly corresponding to 27rn cerl t er , where n center is the 
density at the center of the trap. This may indicate that 
in the trapped system, the CDW order in PSF and SF 
states has a wave vector corresponding to the density at 
the center of the trap. For the CFSF state, because the 
system has a plateau at half filling, the wave vector 2kp 
is 7T /cll. Compared to the homogeneous case, this feature 
is slightly diminished due to the effect of the coexisting 
SF state in the trapped system. In Fig. [8] we show one 
case where CDW order coexists with PSF order in a trap. 
The system size is 80 sites and there are 20 particles of 
each species. The trap frequency £1 = 10~ 5 [/ , t = 0.Q2U 
and U12 = — 0. lit/. The density at the center of the trap 
is roughly 0.45 per site. The cusps are developed around 
±0.97r, which is roughly 2im center . 
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C. Determination of Luttinger parameters from 
experimental data 

Another important question is whether we can use the 
noise correlation Q12 to measure the Luttinger parame- 
ters, K s and K A , in the PSF and CFSF regimes. The 
LL calculation shows that as the system size approaches 
infinity, the noise correlation £/ 12 approaches a power law 
decay with the power —1/ Ks in the PSF regime and with 
— 1/Ka in the CFSF regime (see Eqs. 21 and 22 1. In our 
numerical results for Ci2(q) and D 12(g), we indeed find 
that the decay from the peak at q = satisfies the alge- 
braic decay. To find out the power of the algebraic decay, 
we fit the function Ci2(<?) in the PSF regime and Di2{q) 
in the CFSF regime with the fitting function, 



F(q) = A\sm(2q)\ 



-l/K 



(30) 



where B is the minimum value of 612(9) or D 12(g) and 
A and K are the fitting parameters. In the PSF case 
(U12/U = -0.11), we find that K is 1.3 ± 0.1. This is 
indeed very close to the value of Ks , which is estimated at 
1.4±0.1 obtained by the algebraic fit of Rs- In the CFSF 
case (U12/U = 0.2), we find that K is roughly 1.48 ±0.1, 
while the value of Ka extracted from the algebraic fit 
of Ra is also at 1.48 ± 0.12. Because of the singularity 
at q = 0, a reasonable values of K can be obtained by 
a simple algebraic decay function, Aq~ x l K + B, around 
small q. This shows that even in a trapped system, one 
can still assume a algebraical relationship predicted in 
the LL theory (Eqs. 21 and 22 1 and estimate the values 
of the Luttinger parameters by studying the power of the 
decay from the peak at q = 0. 



V. CONCLUSIONS 

We have studied the behavior of noise correlations for 
a binary bosonic mixture in optical lattices. We con- 
sider different regions of the phase diagram and we show 
that the noise correlations have different signatures for 
different phases. In particular, we discuss the measure- 
ment of the noise correlations as a means for detection of 
the paired superfluid (PSF) and counter-flow superfluid 
(CFSF) order. Our study of a harmonically trapped sys- 
tem shows that the inhomogeneity modifies the noise cor- 
relation, due to the coexistence of different orders within 
the trap. These modifications can be understood in terms 
of the results for the homogeneous system. What we 
find very encouraging is that even with the presence of 
a trap, the noise correlations still have distinctive fea- 
tures for each order, and the peak structure of the noise 
correlation Q12 in the PSF/CFSF regime still obeys the 
algebraic decay relationship predicted by the LL theory. 
This means that one can use the noise correlation to es- 
timate the Luttinger parameters in these two regimes. 
All these results would be useful for experiments aimed 
at detecting the pairing and anti-pairing orders that can 
exist in ultracold atom systems. 
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